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For one-component volatile fluids governed by dispersion forces an effective interface Hamiltonian, 
derived from a microscopic density functional theory, is used to study complete wetting of geomet- 
rically structured substrates. Also the long range of substrate potentials is explicitly taken into 
account. Four types of geometrical patterns are considered: (i) one-dimensional periodic arrays of 
rectangular or parabolic grooves and (ii) two-dimensional lattices of cylindrical or parabolic pits. We 
present numerical evidence that at the centers of the cavity regions the thicknesses of the adsorbed 
films obey precisely the same geometrical covariance relation, which has been recently reported for 
complete cone and wedge filling. However, this covariance does not hold for the laterally averaged 
wetting film thicknesses. For sufficiently deep cavities with vertical walls and close to liquid-gas 
phase coexistence in the bulk, the film thicknesses exhibit an effective planar scaling regime, which 
as function of undersaturation is characterized by a power law with the common critical exponent 
— 1/3 as for a flat substrate, but with the amplitude depending on the geometrical features. 

PACS numbers: 68.08.Bc, 68.08.-p, 05.70.Np 



I. INTRODUCTION 

It is well known that a non-planar topography of a 
substrate modifies its wetting properties significantly. 
Accordingly, experimental studies of complete wetting 
on substrates patterned by a two-dimensional lattice of 
nanopits a one-dimensional array of wedges with fi- 
nite depths [2[ , and arrays of microscopic nonlinear cusps 
or of semicircular channels [H demonstrate the strong in- 
fluence of nanocavities on the adsorption behavior rel- 
ative to that on flat substrates. Theoretical studies of 
adsorption in infinitely deep fH generalized wedges [j| [f| 
predict that such substrates modify the wetting expo- 
nents describing the divergence of the wetting film thick- 
ness upon approaching liquid-gas coexistence. Although 
nanopatterning of surfaces may result in drastic changes 
of their wettability, recent theoretical studies have re- 
vealed surprising hidden symmetries, or so-called covari- 
ances, which relate local adsorption properties for various 
different substrate geometries. 

An example of such a geometrical covariance relating 
wedge and cone complete filling has been reported in 
Ref. p}. It has been shown that the equilibrium mid- 
point interfacial heights in a cone and a iwedge obey 
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the relation li°\Afj,,a) = lw' (A/i/2, a), where a is the 
substrate tilt angle and the undersaturation A/x > is 
the chemical potential deviation from bulk liquid-gas co- 
existence. This relation is valid for the leading asymp- 
totic behaviors of li% in the limit A/i — * + . 

In Ref. Q we have demonstrated that complete wet- 
ting of substrates patterned by periodic arrays of grooves 
or quadratic lattices of pits (see Fig.[T]) exhibits a geomet- 
rical covariance similar to the one described in Ref. 
Moreover, we have also identified a range of undersatu- 
rations, within which the midpoint interfacial height can 
be described by a single scaling function for all four types 
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FIG. 1: (a), (b) Top views of the considered substrate geome- 
tries, (a) Quadratic lattice with lattice spacing P of identical, 
cylindrical or parabolic pits of radius R. (b) Periodic array of 
rectangular or parabolic grooves of width 2R. (c), (d) Cross 
sections along the plane y = of arrays of rectangular grooves 
or cylindrical pits (c), and parabolic grooves or pits (d). All 
cavities have a finite depth D. The grooves have a macro- 
scopic extension in the j/-direction. 



of substrate patterns. In this paper we extend the anal- 
yses of Ref. Q, showing inter alia that, in addition to 
the midpoint interfacial height, other parts of the inter- 
facial profile exhibit the aforementioned covariance rela- 
tion, too, while the laterally averaged interfacial profile 
does not. 

Our analysis is based on an effective interface Hamil- 
tonian approach [l^] • This approach has already been 
applied earlier to study wetting phenomena in a sin- 
gle groove [Til, of a periodic array of grooves fl2| . in 
a wedge [ 4 lial . and also of chemically structured sub- 
strates [l5[ • Although this approach does not model 
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the intrinsic structures of the emerging solid-liquid and 
of the depinning liquid-gas interfaces, it provides reliable 
results [10] concerning the growth and the morphology of 
wetting films, phase transitions, etc. It has been shown 
that the leading asymptotic behavior of the effective in- 
terface potential does not depend on the intrinsic struc- 
tures of the emerging solid-liquid and liquid-gas inter- 
faces and is given by the so-called sharp- kink approxima- 
tion for the density distribution (see Sec. [TTJ) [lOj . There- 
fore, this effective interface model is a reliable approach 
for complete wetting considered here, which probes the 
leading asymptotic behavior of the effective interface po- 
tential. 

We consider rectangular or parabolic grooves and 
cylindrical or parabolic pits, taking into account the long 
range of both the substrate potential and the fluid-fluid 
interaction. The temperature T is chosen to be above 
the wetting transition temperature T w of the unstruc- 
tured planar substrate so that for Ap — > + complete 
wetting occurs. (Similar considerations hold for binary 
liquid mixtures upon approaching their demixing tran- 
sition from the mixed phase.) As function of under- 
saturation complete wetting of such substrates reveals 
four different scaling regimes 8]: filling, postfilling, ef- 
fective planar, and planar, with the three crossover val- 
ues between these four neighboring regimes denoted as 
A/Xjy > Ap% > A/ij,. > 0. Below we briefly describe 
all of them, providing also a short summary of the main 
results obtained. The aforementioned covariance relates 
the behavior of the midpoint wetting film thicknesses for 
all geometries and holds within an undersaturation range 
Ap% < A/i < AfJjf^R) (for A/t£ see below, the super- 
scripts p, g refer to pits and grooves, respectively) which 
we call the postfilling scaling regime. In the case of cylin- 
drical pits or rectangular grooves, for A/i \ ApPf^R) 
and for sufficiently large D/R, the analogue of capil- 
lary condensation occurs such that the pits or grooves 
are rapidly, but continuously, filled by the liquid. How- 
ever, in the case of the parabolic pits and grooves, Apfl'^ 
marks the crossover from the continuous power-law fill- 
ing for A/i > A/Ky^ to the postfilling scaling regime. We 
find numerically for R/a > 50 (see Fig. [T]), where a is a 
molecular length scale, Apfff^R) ~ R^ 1 ^ 5 with a small 
positive effective exponent 5, and A/t^ = 2A/iy i; . 

If we denote the equilibrium interface height at the po- 
sition of the symmetry axes of the cylindrical or parabolic 
pits as lp, and in the middle of the rectangular or 

parabolic grooves as l g , we obtain in the postfilling 
regime 

«w>(A M , R) Pi D) = RApg {^L (^j ( i) 

and 

A p (v) = A g (v/2), (2) 

where £/ is a molecular energy scale. Equation @ ex- 
presses the same covariance relation as the one reported 



in Ref. Q. The dimensionless scaling functions A Pj9 (x) 
of the postfilling regime do not depend on D or P, i.e., 
within this regime the midpoint interfacial heights in- 
crease upon decreasing undersaturation in the same way 
for an isolated cavity as for arrays of them. 

For A/i < A/i£ the cavities are completely filled by 
the liquid and the equilibrium local interface heights 
lp,g(x,y) = Ip.g become de facto independent of the lat- 
eral coordinates x = (x, y). In the case of the rectangular 
grooves or cylindrical pits, A/i^j. marks the crossover from 
the postfilling scaling regime to the effective planar scal- 
ing regime within which the wetting behavior of geomet- 
rically patterned substrates can be mapped onto that of 
layered, laterally homogeneous solids. The upper layer of 
those ersatz solids has a thickness D and its composition 
is related to the geometrical parameters R and P. This 
results in the following scaling relations for the effective 
planar scaling regime: 

l p , g (A»,R,P,D) = ($ p , 9 )^(A/i), 

<f P -l--(!) 2 , (3) 

where l w is the thickness of the adsorbed liquid film on the 
corresponding original non-structured planar substrate. 
$ p and <&g are the areal fractions of solid in the top layer, 
z = 0, of substrates with cylindrical pit and rectangular 
groove patterns, respectively. For D — ► the width of 
the effective planar regime, i.e., the range of applicability 
of Eq. [3] vanishes. 

Finally, at A/i = A/t^ ~ D~ 3 with A/t^ < A/t£ the 
systems cross over to the planar scaling regime, in which 
the geometrical patterns are irrelevant. In the case of the 
parabolic pits and grooves, we do not observe the effec- 
tive planar scaling regime. There is rather an extended 
crossover region from the postfilling scaling regime to the 
planar one. 

In Sect.[Tl]we provide connections between microscopic 
density functional theory and the effective interface po- 
tential approach. In Sect. IIIII we describe in detail the 
four scaling regimes, which are revealed by the mid- 
point interfacial heights as functions of undersaturation. 
Sect. IIVI provides a comparison with experimental data 
and Sect. |V] summarizes our main results and provides 
an outlook. 



II. EFFECTIVE INTERFACE HAMILTONIAN 

We start from the following grand canonical density 
functional for inhomogeneous fluids characterized by the 
number density p(r) [161 ]: 

fi[{/>(r)},T,/i] = J d 3 rf HS (p(r),T) 

+ 2 j ^ r J d3r ' w (\ r ~ r 'IM r M r ') 
+ f d\{p s V{v)-p)p{v). (4) 
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The integrations are performed over the region accessi- 
ble to the fluid; w(r) describes the long-ranged attractive 
part of the fluid-fluid interaction potential and decays 
at large distances ~ r~ 6 . The short-ranged repulsion 
is treated according to the Weeks- Chandler- Andersen 
(WCA) approximation [l7l | by introducing a reference 
system of hard spheres with bulk free energy density 
Ihs(p,T); p is the chemical potential and V(r) is the 
substrate potential, which can be obtained by integrating 
the fluid- substrate pair potential w s (r) over the region 
occupied by the substrate particles with number density 
p s . The equilibrium density distribution minimizes f2. In 
order to make analytic progress we seek the minimum of 
n[p(r)] within the subspace of steplike varying density 
profiles HHol: 

p(r=(x,z)) = 6fz-s(x= (x,y)n 

x ( Pl e(l(x)-z)+p ff e(«-Z(x))),(5) 

where pi and p g are the number densities of the bulk liq- 
uid and gas phases, respectively; s(x) denotes the sub- 
strate surface which is in contact with the fluid; z = Z(x) 
is the local position of the liquid-gas interface assuming 
the absence of bubbles and overhangs; and 6 is the Heav- 
iside step function. Inserting Eq. ((SJ) into Eq. Q yields 
the grand canonical potential f2 as a functional of the 
interface morphology l(x): 

fi[{Z(x)}; pi, p g , fjb] = const 

d 2 z(AQZ(x) + E Iff (x,Z(x)) + W(x,Z(x))). (6) 

A is a macroscopic area, const denotes the sum of all 
terms independent of /(x), and AQ is the difference of 
the grand canonical potential densities of the uniform 
bulk liquid and gas phases, i.e., Ail = Qt,(pi,T, //) — 
D,b(p g ,T, //). Close to the bulk liquid-gas coexistence 
line p = po(T) one has Aft — ApAp + 0{Ap 2 ), with 
Ap = po(T) - p, < 0, and Ap = pi - p g . £; 9 (x, Z(x)) 
describes the cost in free energy to maintain a laterally 
inhomogeneous interface configuration relative to that for 
the flat reference configuration (see Ref. [HI): 

£ iff (x,Z(x)) = 

oo i(x)-/(x') 

-(Ap) 2 J d 2 x' J dz J tfe'tuflr - r'|). (7) 

AO 

The functional form of the effective interface potential 
W(x, Z(x)) depends on the substrate shape as well as the 
fluid-fluid and the substrate-fluid interaction potentials: 

oo oo 

W(x,l(x)) = Ap J d 2 x' J dz J dz 1 

A Z(x)-s(x') 2 

x (piw{y/ (x - x') 2 + z' 2 ) - p s w s {^ (x - x') 2 + z' 2 )) .(8) 
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FIG. 2: Liquid-gas bulk phase diagram in terms of T* = 
fesT/e/ and p* = pa 3 . The critical point is given by 
CT*,p;)« (1.4124,0.2473). 



The fluid particles are assumed to interact via a Lennard- 
Jones potential 



(r) = 4s } 



(- 

V r 



12 /a\' 
r J \r 



(9) 



We apply the WCA [l7| procedure to split up the inter- 
action into a repulsive part 4> rep (r) and an attractive part 
<patt(r), where 



and 



<t>re v { r ) 



4>att{r) 



(r)+s f , r<2 1 ' 6 a 
r > 2 1 /V 



>(r). 



r < 2 1 / 6 cr 
r > 2 1 / 6 cr. 



(10) 



(11) 



The repulsive part gives rise to an effective, temperature 
dependent, hard sphere diameter 



2 1 / < V 



d(T) = J dr 1 - exp( 



4> rep (r) 



(12) 



This expression is inserted into the Carnahan-Starling 
approximation for the free energy density of the hard- 
sphere fluid [HI 

f HS ( P ,T) = k B Tp(\n(p\ a ) + (13) 

r, = \-K P d{T)\ 
6 

where A is the thermal de Broglie wavelength. The at- 
tractive part 0Qtt(r) of the fluid-fluid interaction poten- 
tial is approximated by 



4cr 3 



w{r) = w Q ^-(cr 2 +r 2 ) 3 , 



(14) 
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such that the integrated strength of w(r) equals that of 
the original <j> att {r): 

w Q = fd 3 rw(r) = [ d?rcp att {r)=-^^e f a\ (15) 



Equation (Tl4[) resembles, but differs from the WCA ex- 
pression <j> a tt(f) for the attractive potential; this approx- 
imation offers the benefit to be able to perform certain 
integrations analytically. The actual values pi and p g 
minimize flt(p, T, fj,), and the coexistence line p — po(T) 
is determined by the following equations: 



dp 



dp 







P=Ps 



and 



flb(pi) = Qb(Pg)- 



(16) 



(17) 



The grand canonical potential density of the bulk system 
has the form 

n b (p, T, p) = f HS (p, T) + ^w p 2 - pp. (18) 

The bulk phase diagram resulting from Eqs. (flol - fTS)) is 
shown in Fig. [21 

Since here we are only interested in complete wetting 
involving sufficiently thick films, the functional form in 
Eq. (TT4")) can also be used to model the interaction poten- 
tial between fluid and substrate particles, because under 
these circumstances only its asymptotic behavior at large 
distances matters: 



128\/2 
9tt 



2\-3 



(19) 



In the following we use a s = a. This leads to the follow- 
ing expression for the effective interface potential: 

W(x,I(x)) = ix/ y ,(x,!(x)), 

OO 

Jy s (x,Z(x)) = J dzj dV(a 2 + |r-r'| 2 )- 3 , (20) 

i(x) V, 

where A = - 12 ^ a 6 Ap(p l ef - p s e s ) > plays the role 
of an effective Hamaker constant, and V s denotes the 
domain occupied by the substrate particles which are as- 
sumed to be homogeneously distributed within V s with 
number density p s . The effective interface potential of 
the planar substrate occupying the region z < is 



TTr .,. Air / irl , , / . 

W n U) = — t [a — — + I arctanf — ) 
4a- 3 V 2 a 



(21) 



with W„(l > a) w Air/{121 2 ) so that l v {Ap -> 0) = 
(Air I (GApAp)) 1 / 3 . By applying a small gradient ex- 
pansion [18| for the nonlocal expression E/ g (x, Z(x)) one 



obtains the familiar local approximation for the grand 
canonical potential (see Eq. ©): 

fl loc [{l{x)};pi,p g ,p} = const 
+ J d 2 x(oi g ^J\ + {Vl) 2 + ApApl + W(*,l)y,{22) 

o~ig is the surface tension of the free liquid- vapor interface, 
which within the present sharp-kink approximation has 
the following form: 



OO OO 



(Tig 



-~(A/?) 2 J dz J dz' J d 2 w(\/x 2 + z' 2 ) 



z A 

mV2{Apa 3 ) 2 e f 
9^2 • 



(23) 



Having determined the bulk coexistence curve /j,q(T), 
in the following we minimize the functional in Eq. (|22p 
numerically which yields the equilibrium interface height 
Z(x) within mean-field approximation. The bulk liq- 
uid and gas densities pi and p g , respectively, are eval- 
uated at coexistence, i.e., pi = pi(T, po(T)) and p g — 
p g (T, p (T)). All results presented below are obtained for 
T/T c = 0.85. For the substrate structured with grooves 
the four integrals in the Eq. ([2D]) can be calculated ana- 
lytically, while for the two-dimensional lattices of cylin- 
drical or parabolic pits only three of them can be carried 
out analytically. Whereas the capillary wave spectrum 
of Q given by Eq. ([6]) differs qualitatively [3, [2(| from 
that of the local approximation in Eq. (f2"2"| . the morphol- 
ogy of the wetting film is captured reliably by Eq. (|2"2"]) 
apart from regions of high curvature [2~J ] which do not 
play an important role in the present context. Therefore, 
it is legitimate to apply Eq. (|2"2"]l even in the presence of 
long-ranged fluid-fluid interactions, which are addition- 
ally encoded in W(x, I). 



III. NUMERICAL RESULTS 

For grooves, in Fig. [3] typical interfacial profiles l g (x) 
are shown as obtained for several values of Ap. We are 
primarily concerned with the behavior of the midpoint 
heights as functions of the undersaturation Ap. For 

all substrate geometries the midpoint heights l^l(Ap) 
exhibit four different regimes. 



A. Large Ap: filling and postfllling regimes 

The first regime corresponds to the filling of the cav- 
ities. Cavities with vertical walls (rectangular grooves 
or cylindrical pits) are filled in a qualitatively different 
way compared to parabolic ones. For the case of cylin- 
drical pits and rectangular grooves and for D/R large 
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FIG. 3: Interfacial profiles l g (x) (a) for rectangular grooves 
(which are translationally invariant in the y direction) and 
four values of undersaturation (top to bottom): Afi/ef = 
0.00921, 0.00928, 0.00929, 0.0939 which are close to Afj, 9 fil , 
R/a = 123; (b) for parabolic grooves for Afi/e f = 0.008, 0.01, 
0.012, 0.02 (top to bottom), R/a = 120. For both structures 
P/a = 250 and D/a = 500 (see Fig. []}. I = corresponds to 
the top of the substrate structure at z = 0. Here the effective 
Hamaker constant is A/sf ~ 0.330, and the surface tension 
of the free liquid-gas interface is ai g a 2 /ef ~ 0.478. 



FIG. 4: (a) Interfacial heights (measured from the top sub- 
strate layer at z = 0) at the middle of the cylindrical pits, 
ip ' (dashed lines), and of the rectangular grooves, 4 0) (solid 
lines), as a function of undersaturation A/t measured in units 
of the energy scale £/ of the fluid-fluid interaction, (b) Rescal- 
ing the variables according to Eqs. ((TJ) and ((2| leads to data 
collapse in the postfilling scaling regime; C p = 1 and C g = 2. 
The depth of the cavities is D = 500. All length are mea- 
sured in units of a. A/sf ~ 0.330, and ai g a 2 /ef ~ 0.478. 
A(x) denotes the scaling function introduced in Eq. ([2}. 



enough, there exists a well denned value of the undersat- 
uration A/i = Aji^ at which an abrupt, but still contin- 
uous, filling of the cavities takes place. A similar behav- 
ior has been reported earlier for an isolated rectangular 
groove [ll[. The filling of the vertical cavities is shown 
in Fig. IUa), where Zp°g(A/i) are shown for several values 
of R and P. The interfacial profiles, which are shown 
in Fig. [3(a), are calculated for undersaturations close to 
Afj, 9 fil , | A/i - A/rU < 0.002, and R/a = 123. de- 
pend only on the lateral size R of the cavities, but not on 
the depth D. We find numerically that A/Af ~ i? _1_<5 , 
with an effective exponent 6 ~ 0.035 for both the cylin- 
drical pits and the rectangular grooves. 

Different from the previous case, the complete fill- 
ing of the parabolic cavities is described by an effective 
power law, 4°s(A/i, R, D) ~ A/i K- R > D ) j which is valid 
for A/i > A/tf'y only (see Fig. [5(a)). The effective ex- 
ponent 7 ranges from ca. 3.1, for R — 245er, to ca. 2.0, 
for R = 50cr, at a cavity depth D = 500o\ The value 
7 = 2 is expected to hold for an infinitely deep, single 



parabolic wedge [5j . For both vertical and parabolic cav- 
ities, the complete filling obeys the covariance relation 
4 0) (A/t, R, D) = 4 0) (A/i/2, R, D) valid for A/t > A/x^f. 

For A/i < Afjjft the midpoint height for all patterns 
shows an almost linear dependence on A/i on normal 
scales. The interfacial profiles still reflect the geomet- 
rical patterns, i.e., there are considerable lateral varia- 
tions of the interfacial heights. For large enough D there 
exists a range of undersaturations (Afi^., AfXjf^ within 
which the midpoint interfacial height for all four sub- 
strate geometries can be described by a single scaling 
function We call this range of A/i the postfilling 
scaling regime. The lower bound A/i^ of this regime is a 
decreasing function of the cavities size R. We find that 

for A/i 6 (A/i£, A/t^) the slopes of the 4>°g curves scale 
as R a with a w 2. Combining this result with the scaling 
of the chemical potential A/i^j at the crossover to the 

filling scaling regime, we propose for the functions l^l in 
the postfilling regime the scaling forms given by Eq. fT]). 
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FIG. 5: (a) Interfacial heights at the middle of the parabolic 
pit, 4°' (dashed lines), and of the parabolic grooves, l g ^ 
(solid lines), as a function of A/i. (b) Rescaling the vari- 
ables according to Eqs. ((TJ and ((2} leads to data collapse for 



A M G (Aa4,A^) 



C 71 



1 and C 9 = 2. The depth of the 



cavities is D — 500. All length are measured in units of a. 
A/ej ~ 0.330, and ai g a 2 /£f » 0.478. A(a;) denotes the scal- 
ing function introduced in Eq. ([2]). The filling regime ends at 
large A/i (dotted line) where the wetting films in the center 
become microscopicly thin so that scaling breaks down. 



In this postfilling scaling regime the scaling function A g 
describes rectangular as well as parabolic grooves, and 
similarly A p describes cylindrical and parabolic pits. 
The scaling functions A p s and the corresponding data 
collapse upon suitably rescaling the chemical potential 
and the midpoint heights are shown in Figs. HJb) and 
5(b) for the case of vertical and parabolic cavities, respec- 
tively. In accordance with Eq. the difference between 
A p (x) and A ff (a;) can be absorbed in the rescaling factors 
C p = 1 and C g = 2. Inspection tells that the resulting 
common scaling function denoted as A(x) in Figs. HJb) 
and[5Jb) are indeed identical. The scaling functions A Pj9 
in the postfilling scaling regime do not depend on the 
cavity depth D which is demonstrated in Fig. [6l where 
Zp°^(A/i) is shown in the case of cylindrical pits and for 
several values of D. For even smaller values of D/R, ow- 
ing to the finite-size rounding effects, the filling of the 
cavities does not occur at a certain value of the chemical 
potential, but there is an entire transition region. As a 
result, the range of the postfilling scaling regime shrinks 



FIG. 6: Interfacial heights 4°' m the middle of the cylindrical 
pits as a function of undersaturation A/i for several values 
of the pit depth D (see Fig. [TJ. All lengths are measured in 
units of a; A/ef w 0.330 and ai g a 2 /e f w 0.478. 




FIG. 7: Interfacial profiles within the cavity regions for 
cylindrical pits, l p (x,y = 0) (dashed lines), and rectangular 
grooves, l g (x) (full lines), for the postfilling scaling regime. 
Black, red, and blue colors correspond to (R — 50, P = 110), 
(R = 120, P = 250), and (R = 245, P = 500), respectively. 
The dashed and full interfacial profiles are calculated for two 
different values of the undersaturation which both satisfy 
C p ,gAn/s f (R/a) im5 = 1.18 with C p = 1 and C g = 2, re- 
spectively. In the inset the axes are rescaled by _R _1 leading 
to complete data collapse. D — 500 and a is the length unit; 
A/ef » 0.330 and a lg a 2 /e f « 0.478. 



as D/R decreases. A PiS do not depend on the pattern 
periodicity P either. In the case of a single cavity, i.e., in 
the limit P — * oo, we obtain the same curves l p °l(A^) as 
those presented in Figs. H] and [5] Furthermore, we find 
that for the filling and postfilling scaling regimes not only 

l p °l, but also the full interfacial profiles i P)5 (x) are, within 
numerical accuracy, independent of the pattern period- 
icities P studied here. This implies that for the case of 
pits (both cylindrical and parabolic) the interfacial pro- 
file within the cavity region is, to a good approximation, 
rotational symmetric around the symmetry axis of the 
pit, i.e., Z p (x) = l p (r p ), where r p < R is the radial dis- 



tance from the symmetry axis of a pit. 

In the postfilling scaling regime, in addition to 



the 
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midpoint interfacial heights, those values of the profiles 
l Pt g(x), which are taken off center at the same distance 
r Pt g/R < 1 from the symmetry axis of a pit or from the 
midplane of a groove, also satisfy generalized scaling re- 
lations similar to Eqs. (P) and ([2]): 

l„(x, A/i, R, P, D) = RA p (^, ^ (?) 1+5 ^j (24) 

l g (x, Afi, R, P, D) = RA g (!±, ^(?) 1+S ), (25) 

where r p — y/(x — Xq) 2 + (y — j/o) 2 is the radial distance 
from the symmetry axis of a pit at (x,y) = (xq, j/o); r g = 
|x — xo| is the distance from the midplane of a groove 
at x — x . The scaling functions A p and A g obey the 
covariance relation 

A p (u,v) = A g (u,v/2), u<l. (26) 

The interfacial profiles for cylindrical pits, l p , and rect- 
angular grooves, l gi are shown in the Fig. [7] for several 
R. The profiles are calculate for two different under- 
saturations Ap within the postfilling regime which both 
satisfy C P}g Ap/ e } {R/ a) lmb = 1.18 with C p = 1 and 
C g = 2, respectively. Rescaling the variables according 
to Eq. ([25)) leads to data collapse. Equation (|26|) means 
that, for fixed R and certain values of undersaturation, 
the interfacial profile within a pit region can be consid- 
ered as a surface of revolution generated by the function 
l g (r g ) which characterizes the interfacial profile within a 
groove region. This also holds within the filling scaling 
regime. 

We also consider the lateral average l p . g = 
J d 2 x l Pi g (x) I A of the interfacial profiles l Pi9 (x) with the 
integral running over the lateral surface area A of a unit 
cell of the surface pattern. This quantity is related to 
the excess adsorption r = J d 3 r(p(r) — Pb)/A per area of 
the cell, with the integral running over the region acces- 
sible to the fluid above a unit cell. Using the sharp-kink 
approximation for the density (see Eq. T can be ex- 
pressed as 

T p , g = jAp J ^(/(x) - s(x)) 
A 

= + (27) 

where V p ^ g is the volume of a pit or groove, respectively. 
For rectangular grooves and cylindrical pits Fig. [8j shows 
the corresponding results for / Pj£) (A/i) with the chemical 
potential and I rescaled according to Eqs. U) and (J5J. 
The laterally averaged interfacial profile in the postfill- 
ing scaling regime also obeys the scaling described by 
Eq. |T]), but the scaling functions for the rectangular 
grooves and the cylindrical pits do not satisfy the covari- 
ance relation given by Eq. ([2|) . The same result holds for 
the case of parabolic pits and grooves, which is shown in 
Fig. Moreover, the scaling functions which are shown 
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FIG. 8: Rescaled laterally averaged interfacial profile for 
cylindrical pits, l v (dashed lines) , and for rectangular grooves, 
l g (full lines) , as a function of undersaturation Afi. The rescal- 
ing of the variables according to Eqs. |T} and ([2]) leads to two 
different scaling curves describing the cylindrical pits and the 
rectangular grooves, respectively. Thus the covariance rela- 
tion given by Eq. is not fulfilled. Compare Fig. |3£b) which 
shows that the postfilling scaling functions for the midpoint 
interfacial height are the same for these two cavity types. 
C p — 1 and C g = 2, D — 500, and a is the length unit; 
A/e f w 0.330 and ai g a 2 /£ f w 0.478. 



in Fig. [51 describe also the laterally averaged interfacial 
profiles for parabolic pits and grooves, respectively. Ge- 
ometrical considerations show that the lateral average of 
a troughlike surface generated by a function of a single 
variable is smaller than the lateral average of the cor- 
responding rotationally symmetric surface generated by 
revolving the same function around the symmetry axis. 
As shown above, both surfaces obey the covariance scal- 
ing properties. Therefore this inequality of the lateral 
averages implies that this covariance does not hold for 
the lateral averages. 



B. Small A/j, : effective planar and planar scaling 
regimes 

For undersaturations below a certain value A/^j., i.e., 
closer to gas-liquid coexistence, the liquid-gas interface 
becomes flat (Fig. [TU]) . its height denoted as l p . g for the 
pit and groove patterns, respectively. For these flat wet- 
ting films the functional in Eq. (f2"2"| reduces to a function 
of a single variable l Pjg the minimum of which determines 
the equilibrium thicknesses of the wetting films in this 
regime. The corresponding results l p . g (Ap) for cylindri- 
cal pits and rectangular grooves are presented in Fig. [11] 
and Fig.[T21 respectively. For these two substrate geome- 
tries there is an interval of undersaturations (A/i^, A//£), 
within which both l p and l g reveal the same power law 
growth l Pi g ~ A/i -1 / 3 as for a flat surface, but with differ- 
ent amplitudes reflecting different, geometry dependent, 
effective Hamaker constants A e . We call this range of un- 
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FIG. 9: Rescaled laterally averaged interfacial profile for 
parabolic pits, l p (dashed lines), and grooves, l g (full lines), 
as a function of undersaturation Ap. The rescaling of the 
variables according to Eqs. ([T]) and © leads to two different 
scaling curves describing the cylindrical pits and the rectan- 
gular grooves, respectively. Compare Fig. [5jb) which shows 
that the postfilling scaling functions for the midpoint interfa- 
cial height are the same for these two cavity types. C p = 1 
and C g = 2, D = 500, and a is the length unit; A/sf « 0.330 
and ai g a 2 /e f w 0.478. 




(a) (b) 



FIG. 10: (a) Schematically drawn liquid wetting film mor- 
phology on substrates sculptured by cylindrical pits or rect- 
angular grooves within the effective planar scaling regime, (b) 
Wetting film of thickness 1% on a flat layered solid. The top 
layer has a thickness D with an effective Hamaker constant 
A e whereas the semi-infinite bottom part of the substrate is 
characterized by the Hamaker constant A. 

dersaturation the effective planar scaling regime. Within 
this regime the substrates sculptured by cavities with ver- 
tical walls are equivalent to layered and flat ersatz solids 
(see Fig. [TTiT b)) with the effective interface potential fl2| 

it { A A e A e \ 

Here the first term corresponds to a flat semi-infinite, 
z G (— oo, —D), solid with Hamaker constant A and the 
second and the third are the corrections due to the finite 
slab, z G (—D,0), with effective Hamaker constant A e . 
For a <C I < D, it is a good approximation to ignore the 
bottom part of the substrate (i.e., the first and the third 
terms in W(l 3> a)) so that the amplitude is determined 
by the effective Hamaker constant A e . 

With A s we denote that region of the sculptured part 



FIG. 11: Wetting film thickness l p above cylindrical pits 
for small Au. The solid lines are obtained by minimiz- 
ing Eq. (1221) in the subspace of laterally constant interfacial 
heights Z(x) = l p . The circles represent the mid-point in- 
terfacial heights ip°' and are the results of the full numerical 
minimization of the Hamiltonian in Eq. (|22p . The circles are 
plotted in that Afi range for which the deviations from the 
full and dashed lines become detectable. The dashed lines 
correspond to the film thicknesses l%(Au) on layered flat er- 
satz substrates (see Fig. 1 1 Of) . with a top layer of height D and 
an effective Hamaker constant A e = A& p = A(l — tt(R/P) 2 ), 
where A is the Hamaker constant of the solid without pits. 
In the inset the vertical axis is rescaled according to the scal- 
ing relation given in Eq. Q leading to data collapse and - 
after rescaling - agreement with l n within an intermediate 
regime, Afi^ < An < A/z|, the width of which increases as 
A/j,% — const x D~ 3 for D —* oo. The pit depth is D = 500cr; 
A/ef ~ 0.330 and a lg a 2 /e f » 0.478. 

of the substrate (z G [— D,0]), which is occupied by the 
solid, and with Ao the entire slab z G [-D, 0]. The equiv- 
alence of the patterned and the flat layered substrate ac- 
tually means Wa,(x,Z) = W^ o (x,Z). With Eq. one 
obtains A e = AI\ s /Ia - For I a the integrals 7a; 
(i — s, 0) can be approximated as 

oo 

f , f dV ( CT 2 + ||X - X'|| 2 \ , . 

W^/f^tH 7? > (29) 

I Ai 

Actually, I\ a is a functions of x, which varies across a unit 
cell and only its leading term ~ l~ 2 is laterally constant, 
fn Eq. taking into account only the leading term as 
function of I one finds A e sw AS s /So, where S s /So is the 
fraction of area occupied by the solid in the top layer of 
the substrates patterned by cavities with vertical walls. 
Thus, one has A e g sa A$> g for the rectangular grooves 
and Ap w A$ p for cylindrical pits (see Eq. flSJ). The 
thicknesses i^(A^t) of the wetting films on such layered 
ersatz substrates, which effectively correspond to arrays 
of grooves, are almost indistinguishable from the corre- 
sponding actual ones l g . For lattices of cylindrical pits in 
Fig. QT]we present /J(A/i) (dashed lines) on layered sub- 
strates, and find agreement with the actual thicknesses 
Ip(Afi) within the effective planar scaling regime. There- 
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FIG. 12: Wetting film thickness l g above rectangular groove 
patterns for small A[i. The solid lines are obtained by min- 
imizing Eq. (|22p in the subspace of laterally constant inter- 
facial heights Z(x) = l g . The circles represent the mid-point 
interfacial heights lg and are the results of the full numer- 
ical minimization of the Hamiltonian in Eq. (|22|l . The cir- 
cles are plotted in that range for which the deviations 
from the full and dashed lines become detectable. On the 
present scales the data for l g are indistinguishable from the 
lines l%(Afi) obtained for layered flat ersatz substrates with 
a top layer of height D and an effective Hamaker constant 
A e — A$ 9 = A(l — 2R/P), where A is the Hamaker constant 
of the solid without grooves. In the inset the vertical axis 
is rescaled according to the scaling relation given in Eq. © 
leading to data collapse and - after rescaling - agreement 
with 1^ within an intermediate regime, A^ < A/j, < A/j.%, 
the width of which increases as Afi% — const x D~ 3 for 
D -> co. The groove depth is D = 500cr; A/e f « 0.330 
and ai g a 2 /e f w 0.478. 



fore we can conclude that for sufficiently thick wetting 
films Ip^g obey the scaling relation given in Eq. ([3]). In- 
deed, after rescaling the film thicknesses by the geometry- 
dependent factors (^p^g)^ 1 ^ 3 , the curves for l p , g collapse 
and, within the numerical precision, they coincide with 
the curve for l„. This is shown in the insets of Figs. [11] 
and [T2j In the case of parabolic cavities we have found 
numerically that the corresponding curves for Z Pi9 do not 
reveal this effective planar scaling regime, but gradually 
cross over to the planar one, denoted as l„ Q. We em- 
phasize that the effective planar regime exists only for 
patterns with vertical walls. Thus, we also expect this 
regime to exist for lattices of cylindrical nanoposts, but 
not for, e.g., sawlike substrate surfaces or periodic arrays 
of trapezoidal cavities. 

As l Pt3 reach the value ~ D, i.e., at = Afi^ ~ D~ 3 , 
a crossover to the planar scaling regime takes place. In 
the planar scaling regime Afj, < A/i^ the geometrical 
patterns have become irrelevant. Thus the range of ap- 
plicability A/i^ — A/in of Eq. [3] increases with the cavity 
depth D as A^ — const x D~ 3 , whereas for D — > Afi v 
merges with A/i^., eliminating the effective planar scaling 
regime. 



IV. COMPARISON WITH EXPERIMENTAL 
DATA 

In this section we compare our results with the corre- 
sponding experimental data of Ref . |l| , where the adsorp- 
tion of methylcyclohexane (MCH) on a nanopatterned 
silicon substrate has been studied using X-ray reflectiv- 
ity (XR) and grazing incident diffraction (GID) tech- 
niques. The substrate was patterned by a hexagonal 
lattice of parabolic pits with radius R « 123 A , depth 
D w 200A , and lattice constant P « 394A. Accord- 
ingly, in the following the theoretical data have been ob- 
tained for a hexagonal pattern, too. The data in Ref. 1] 
have been interpreted in terms of the excess adsorption 
r = Ti n + T ou t as the sum of the excess adsorption r.;„ 
in the pit and the excess adsorption r oMt above the pit 
opening (see inset in Fig. EUa)). These two quantities 
are defined as 



o 



and 




Using the sharp-kink approximation for the density pro- 
file (see Eq. ©) Eqs. J30]) and ^ take the form 



and 




In these experiments the thermodynamic deviation from 
liquid- vapor coexistence of MCH has been tuned via con- 
trolling a temperature difference AT between the sub- 
strate and the reservoir of bulk liquid. AT is related to 
the undersaturation A/i according to [l| AT = A/iT/H, 
where T = 305K is the reservoir temperature, and 
H = 5.8 x 1CP 20 J is the latent heat of vaporization 
per MCH molecule We have modeled MCH by a 
Lennard-Jones interaction potential with the parameters 
a = 5.511A and e f /k B = U6K [13], which leads to 
AT w 32.37-^X. By fitting the experimental reference 
data for the wetting film thickness on the corresponding 
planar substrate [1] to the expression Z ff (A/x — > 0) given 
just below Eq. (|2"Tj) we are able to determine the ratio 
A/ (e fAp*) « 1.006 involving the Hamaker constant A. 
The last material parameter required to facilitate a quan- 
titative comparison between our theory and the experi- 
mental data in Ref. [l[ is the liquid-vapor surface tension 
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FIG. 13: Excess adsorptions (a) T out (AT) above the top 
of the substrate and (b) Tin (AT) within the parabolic pits. 
Symbols and black lines denote experimental data from 
Ref. [l|]. Red lines are the corresponding theoretical results 
(see main text). p e = 0.26e/A 3 is the electron density of the 
bulk MCH. AT is proportional to A/z. 



of MCH, for which we use oi„ = 22.72^ for T = 30°C 

■ 1 ' L y cm 

|23l ]. Using the values of a and £/ given above for MCH, 
we find a lg a 2 /e f « 1.149. According to Eq. IJ23J) this 
value of the surface tension corresponds to Ap* w 0.668 
or T* w 0.88 as follows from the bulk phase digram shown 
in Fig. [2j It is satisfactory to see that this value of the 
reduced temperature T*/T* = T/T c w 0.88/1.41 w 0.62 
turns out to be reasonably close to the actual experimen- 
tal ratio T/T c = 305/572 » 0.53. With this there are no 
free parameters left in the model. 

T out and Ti n as functions of AT are shown in Fig. ll3f a) 
and Fig. fT3Fb) , respectively. The X-ray reflectivity tech- 
nique, which has been used in Ref. [l[ , provides access to 
the laterally averaged electron-density profile p(z); conse- 
quently the adsorption is measured in units of the areal 
electric charge density ek~ 2 . In order to present the 
theoretical adsorption data in these units, we use for 
the electron-density of the bulk liquid MCH the value 
p e = 0.26eA~ 3 and neglect the density of the vapor. For 
temperature differences in the range 5K < AT < 8K 
we identify the filling regime, i.e., within this range the 
theoretical data I\- n (AT) resemble a power-law growth 
Ti n oc AT~P p with an effective exponent (3 P ~ 2.4. This 



effective exponent (3 P turns out to increase with the pit 
size and for R/a — 245 and D/a — 500 we find (3 P ~ 3.4. 
There is no obvious connection between the effective ex- 
ponent f3 p and the exponent j(R, D), which is introduced 
in Subsec. HIT Al for the midpoint interfacial height. In 
Ref. [l[ the value of effective exponent (3 P ~ 3.4 has been 
roughly estimated by using 7 = 2 (which is valid for in- 
finitely deep parabolic wedges or pits) and the expression 
V(l p 0) ) = l(Rl p 0) ) 2 /D + f (R 2 l p 0) /D) 3 / 2 for the volume 
enclosed between the surface of parabolic pit and a spher- 
ical meniscus 0, [j| which approximates the shape of the 
interfacial profile. With this, for the experimentally rel- 
evant values of l p °^ the dependence of V on A/i can be 
approximated as V(l p ^ ~ Ap~ 2 ) ~ A^ -3 - 4 [l|. 

At AT w 8K one can identify a small change in slope 
of the (AT) curve. This crossover corresponds to the 
disappearance of the interfacial meniscus, such that for 
larger values of AT the interface follows the substrate 
shape. The authors of Ref. [l| assigned the value (3 P ~ 
0.76 to their data. We explain this disagreement with 
the value (3 P w 2.4 reported here (see Fig. [T^b)) by a 
mislocation of the filling regime for the experimental data 
of r in (AT) induced by the large error bars for the data 
points at large AT. There also seems to be a difference 
between the behavior of the theoretical and experimental 
results for T out at large AT (see Fig.lTBTa)). However, for 
such large undersaturations the wetting film is only ca. 
10A thick causing rather large experimental error bars 
(actually larger [2J] than those presented in Fig. ITSTa) 
which are taken from Ref. [l|) for extracting values for 
the adsorption T out from the scattering data. On the 
other hand, for such thin wetting films the reliability of 
our theoretical description is expected to deteriorate, too. 
Due to packing effects the use of an effective interface 
Hamiltonian should be replaced by a full-fledged density 
functional approach. 

In Fig. [T3] the laterally averaged interfacial profile 
l p (AT) is shown. In order to obtain the experimental 
curve l p , first we have interpolated the experimental data 
for Ti n and T out , second we have formed T = Ti„ + T out , 
and by using Eq. (j2"T|) we have determined l p (AT). This 
experimental curve is compared with the corresponding 
theoretical one. The theoretical curve for l p reveals the 
filling and postfilling scaling regimes, as well as the ex- 
tended crossover to the planar scaling regime. Due to 
the relatively small depth D/a ~ 36 of the pits the 
crossover region between the first two regimes is rather 
extended. In the postfilling scaling regime, i.e., for 
0.7K < AT < 4K, we obtain very good agreement be- 
tween the experimental and theoretical data, while there 
are substantial deviations in the filling regime. With re- 
spect to the latter one we suggest that for deeper pits the 
experimental identification of the filling regime would be 
easier and the error bars smaller. First under this con- 
ditions the filling regime occurs over a wider range of 
undersaturations; second, the midpoint interfacial height 
attains larger values; third, the crossover to the postfill- 
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FIG. 14: Laterally averaged interfacial profile l p (AT) for 
parabolic pits, obtained according to Eq. 1271 The black line 
interpolates the experimental data from Ref. [l| (without er- 
ror bars), and the red line denotes corresponding theoretical 
results. AT is proportional to Afj, (see the main text after 
Eq. (133). 

ing regime becomes more pronounced. 

V. SUMMARY AND OUTLOOK 

Within an effective interface Hamiltonian approach 
based on density functional theory (Sect. QT} we have 
studied complete wetting of four classes of solid sub- 
strates topographically patterned on the nanoscale: 
quadratic (and hexagonal) lattices of cylindrical or 
parabolic pits and periodic arrays of rectangular or 
parabolic grooves (Fig. [J). The long-ranged fluid-fluid in- 
teraction (giving rise to the bulk phase diagram shown in 
Fig. and the fluid-substrate potential are modelled by 
isotropic pair potentials decaying with distance ~ 1/r 6 . 
By analyzing the full range of undersaturations A/z rela- 
tive to liquid-vapor coexistence in the bulk we have iden- 
tified four different regimes for the increase of the wetting 
film thickness (Subsecs. IIIII A and IIIII B). and the exis- 
tence of so-called covariances relating this increase for 
pits and grooves (Eq. ©). 

The filling of a cavity (Fig. [3]) by the liquid occurs 
independently of the presence of other ones. This sug- 
gests that, e.g., in the filling regime the power law in- 
crease of the excess coverage in parabolic cavities is not 
only accessible to X-ray scattering (which requires peri- 
odic arrays) but also to AFM techniques applied to single 
cavities (see, e.g., Ref. (H). For D / R £ °- 5 ( see Fig. [J) 
there exists a range of undersaturations, which we call 



the postfilling scaling regime, in which the local interfa- 
cial characteristics, obtained for all four patterns, can be 
expressed in terms of a single scaling function (Eqs. |T]), 
©, ((2SJ), and (HU) and Figs. H EE and[7J). However, the 
characterization of the laterally averaged interfacial pro- 
file requires two different scaling functions, one for pits 
and one for grooves (Figs. [5] and [5]). All these scaling 
functions are independent of the cavity depth (Fig. [5]). 

For small undersaturations, the increase of the wetting 
film thickness is determined by the entire array of cavi- 
ties. In the case of vertical cavity walls, i.e., for rectangu- 
lar grooves and cylindrical pits, the interfacial thickness 
increases as on a planar substrate, but with an effec- 
tive, geometry-dependent, Hamaker constant (Eq. [3] and 
Figs. [TO1TT1I)) . For very small undersaturations the pat- 
terns on the substrate become irrelevant, and the wet- 
ting film thickness increases as if the substrate would be 
flat. The crossover to this latter so-called planar scal- 
ing regime occurs only for long-ranged dispersion forces. 
For short-ranged interactions, instead, the increase of the 
wetting film thickness would remain determined by the 
areal fraction of the solid at the substrate surface z = 
for all film thicknesses [26[ . 

We have been able to report quantitative agreement 
between our theory and the X— ray scattering data of 
Ref. d (Figs. Ql and [SJ). For the range of undersatura- 
tions corresponding to the postfilling scaling regime, the 
agreement is satisfactory (Fig.rjJJ. We have been able to 
clarify the interpretation of the adsorption parabolic pits 
in terms of an effective filling exponent (3 P (Fig. [TWb)). 
We suggest, that for deeper pits the experimental iden- 
tification of the filling regime would be easier. In order 
to test our scaling predictions (Eqs. fl} and @) experi- 
ments for the following surface structures are still miss- 
ing: (i) patterns with cavities of different sizes have to 
be considered; (ii) resolving the full morphology of the 
wetting films, or at least determining the midpoint inter- 
facial heights, are highly desirable. 

A natural extension of the work presented here would 
be to study the temperature dependence of the filling 
regime and of the scaling function for the postfilling scal- 
ing regime. Another perspective is to consider substrates 
patterned by, e.g., cylindrical nanoposts. 
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